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Abstract 

With the "bouncing photon" treatment, the gauge invariance of the 
response function of interferometers to scalar gravitational waves (SGWs) 
has been recently demonstrated in its full frequency dependence in three 
different gauges well known in literature, while in previous works such 
invariance was only proved in the low frequencies approximation. After 
a review of scalar tensor theories of gravity, in this paper the analysis of 
the response function to SGWs is generalized in the full angular and fre- 
quency dependences and directly in the gauge of the local observer, which 
is the gauge of a laboratory environment on Earth. The result is used 
to analyse the cross-correlation between the two LIGO interferometers 
in their advanced configuration for a potential detection of a stochastic 
background of SGWs and to release a lower bound for the integration 
time of such detection. A comparison between the response function in- 
troduced in this paper and previous low-frequency approximated ones is 
also performed. Because the lower bounds result very long, we hope in the 
LISA interferometer and in a further growth in the sensitivity of advanced 
projects. 

PACS numbers: 04.80.Nn, 04.30.Nk, 04.50.4-h 



1 Introduction 

Recently, the data analysis of interferometric gravitational waves (GWs) detec- 
tors has been started (for the current status of GWs interferometers see [1]) and 
the scientific community aims at a first direct detection of GWs in next years. 
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Detectors for GWs will be important for a better knowledge of the Universe 
and also to confirm or ruling out the physical consistency of General Relativity 
or of any other theory of gravitation This is because, in the context 

of Extended Theories of Gravity, some differences between General Relativity 
and the others theories can be pointed out starting by the linearized theory of 
gravity [13]- [HI- In this picture, detectors for GWs are in principle sensitive 
also to a hypothetical scalar component of gravitational radiation, that appears 
in extended theories of gravity like scalar-tensor gravity and high order theories 

m-m- 

With the "bouncing photon" analysis (see figure 1) in [13] the gauge invari- 
ance of the response function of interferometers to SGWs has been demonstrated 
in three different gauges well known in literature and for all the frequencies of 
SGWs, while in previous works this gauge invariance was only proved in the low 
frequencies approximation, see for example |16j . In this paper, after a review 
of scalar tensor theories of gravity, which is due to provide a context to bring 
out the relevance of the results, the analysis of the response function for SGWs 
is generalized in the full angular and frequency dependences and directly in the 
gauge of the local observer, which is the gauge of a laboratory environment 
on Earth. The result is used to analyse the cross-correlation between the two 
LIGO interferometers in their advanced configuration for a potential detection 
of a stochastic background of SGWs and to release a lower bound for the in- 
tegration time of such detection. A comparison between the response function 
introduced in this work and previous low-frequency approximated ones is also 
performed. Because the lower bounds result very long, we hope in the LISA 
interferometer and in a further growth in the sensitivity of advanced projects. 

2 A review of scalar gravitational waves from 
scalar-tensor theories of gravity 

If the gravitational Lagrangian is not linear in the curvature invariants the 
Einstein field equations have an order higher than second [2]-[l9]. For this 
reason, such theories are often called higher-order gravitational theories. In the 
most general case they arise from the action [13] 

J d^xV^lFiR, DR, tfi?, a'^R, 0) - ^^•'(b.^,^,, + Cm], (1) 

where F is an unspecified function of curvature invariants and of a scalar 
field (j) and the symbol □ is the D'Alembertian operator. The term Cm is the 
minimally coupled ordinary matter contribution. 

In scalar-tensor theories of gravity both of the metric tensor f;^ and a fun- 
damental scalar field are involved [TKUdKnillHllTnillllll]. Their action 
can be recovered from cq. ([T]) with the choice [13] 

F{R,<P) = f{(l))R-V{(t)) and e -1 . (2) 
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Considering the choice ^ , the most general action of scalar-tensor theories 
of gravity in four dimensions is given by [13[ I17j 



J d^x^[f{cl>)R + \g^''cl,,,(^.,, - + Cra]. (3) 

Choosing 

^ = ^{^)^47W) W{^) = V{<j,{ip)) (4) 

eq. ^ reads 

j d^xV^i^R - ^gt^-^.^^^.^, - W{^) + £,„], (5) 

which is a generalization of the Brans-Dicke theory [53] . 
By varying the action ([5]) in respect to and to the scalar field (p the field 
equations are obtained: 

(6) 

+^{V;t^u - gt.u'^^) + j^g^^W{Lp), 
with associated a Klein - Gordon equation for the scalar field 



We emphasize that in this paper natural units are used: G = 1, c = 1 and 
h=l. 

In the above equations T^™^ is the ordinary stress-energy tensor of the matter 
and G is a dimensional, strictly positive, constant [13]. The Newton constant is 
replaced by the effective coupling 

Ge// = 

which is, in general, different from G 
the scalar field couphng is 

ip = const. = — -. (9) 

The case in which it is w = const, in eqs. ([6]) and ([7]) is the string-dilaton 
gravity [T31 dH [3Ij ■ If one wants to study gravitational waves, the linearized 
theory in vacuum (T^™'' = 0) with a little perturbation of the background has to 
be analyzed. The background is assumed given by the Minkowskian background 
plus Lp — ipo and (po is also assumed to be a minimum for W [T5j : 

W ~ ^aSp^ ~ aSip (10) 



General Relativity is obtained when 
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Putting 

(11) 

and, to first order in h^^i, and Sip, if one calls Rfj,upa , R^J.v and i? the linearized 
quantity which correspond to R^vpa , Rpu and R, the linearized field equations 
are obtained [T3l: 



where 
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TO 



(12) 



(13) 



2w+3 ■ 

The case in which it is = const, and Vl^ = in eqs. ([6|) and ^ is 
the analogous of the Brans-Dickc theory [23|. For a sake of simplicity, we will 
analyse such a case. 

Then, the linearized field equations become 



Rpu - f r/^^ - df^d^i + ?/^^n$ 



□$ = 0. 



(14) 



Let us put 



h = r/^'' V ^-h- 4^, 
with h = Tj^^h^y^ where the inverse transform is the same 

h = v'^'h^u ^-h- 4e 
By putting the first of eqs. (fTB|) in the first of the field eqs. ([H)) we get 



(15) 



(16) 



aJi^. - a^(5"/i„.) - a,(a"/i„ J + Vf^.a^'id^h^p). (17) 

Now, let us consider the gauge transform (Lorcnz condition) 
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h^h' = h + 2a"e„ (18) 

with the condition Dej^ = d^h^y for the parameter e^. We have 

d^h'^, = 0, (19) 
and, omitting the we can rewrite the field equations hke 

nJi^. = (20) 

□e = 0; (21) 
solutions of eqs. (PO)) are plan waves: 

hf_,^ = A^^Ck-) expiik^Xa) + c.c. (22) 

^ = a{lt ) cxp{ik"xa) + c.c. (23) 

Thus, in eqs. ([20)1 and ([2^ the equation and the solution for the tensorial 
waves exactly like in General Relativity (ref. [21]) have been obtained, while 
eqs. (PT|) and are respectively the equation and the solution for the scalar 
massless mode. 

The solutions 122]) and (|23p take the conditions 

/c kiy^ — 

(24) 

fc^A^, - 0, 

which arises respectively from the field equations and from eq. (|19p . 

The first of eqs. ([24]) shows that perturbations have the speed of light, the 
second gives the transverse effect of the tensorial components. 

Fixed the Lorentz gauge another transformation with De^ = can be per- 
formed; let us take 

De^' = 

(25) 

which is permitted because = = Oh. One gets 

h^^2^ ^ V = V . (26) 

i.e. /i^jy is a transverse plane wave too. One can also see that the gauge 
transformations 
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□e^ = 

(27) 

5^6'' = 0, 

save the conditions 

df'h^, = 

(28) 

h = -2^. 

Consider a wave incoming in the positive z direction it is 

fc^ = (fc,0,Ofc) (29) 
and the second of eqs. ([M)) imphes 

A,o = -A,3 (30) 

^00 = ^^30 + ^33- 

Now let us compute the freedom degrees of A^^. We was started with 
10 components (A^jy is a symmetric tensor); 3 components have been lost for 
transverse effect, more, the condition (j26p reduces the component to 6. One can 
take Aqq, All, A22, A21, A31, A32 like independent components; another gauge 
freedom is present }13j . and it can put to zero three more components (i.e. one 
can only chose three of e^, the fourth component depends from the others by 
a^e'^ = 0). 



Then, taking 



ep( k ) eKp{ik°'Xa) + c.c. 



(31) 



the transform law for A^j,^ is (see eqs. and ([^ ) 



Afj,u > A^j^ — A^j_i, 2ifc(^ejy). (32) 
Thus, for the six components of interest 



^00 Aqo + 2ifceo 

All All 

A22 — * A22 

A21 A21 

A31 A31 - ikei 

A32 A32 - ikt2. 



(33) 



The physical components of A^^, are the gauge- invariants An, A22 and A21, 
thus one can chose to put equal to zero the others. 
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The scalar field is obtained by eq. ([^5]) : 

h^h^hn+ h22 = -2^. (34) 

Now, defining $ = — ^, the total perturbation of a gravitational wave incom- 
ing in the z+ direction in this gauge is 

V(t -z) = A+it- z)eS;j + AX(t - z)e(x; + <f{t - (35) 

The term (t — z)e'lit^ + (t — z)e^^^ describes the two standard (i.e. ten- 
sorial) polarizations of gravitational waves which arises from General Relativity 
in the TT gauge [H], while the term $(t — 2:)ejfj is the extension of the TT 
gauge to the scalar case. 

For a purely scalar GW the metric perturbation ((35l) reduces to 

K. = ^el'l (36) 
and the correspondent line element is |13| 

ds^ = dt^ - dz^ - (1 + $)da;2 - (1 + $)dy^ (37) 

In other words, in scalar-tensor theories of gravity, the scalar field generates 
a third polarization for gravitational waves. In the analyzed case three differ- 
ent freedom degrees are present, differently from the standard case of General 
Relativity, where the freedom degrees for gravitational waves are two [23] 

3 Analysis in the gauge of the local observer 

As the potential detection of GWs is performed in a laboratory environment on 
Earth, the coordinate system in which the space-time is locally flat is typically 
used and the distance between any two points is given simply by the difference 
in their coordinates in the sense of Newtonian physics [131 [24j [25] . In this gauge, 
called the gauge of the local observer, SGWs manifest themselves by exerting 
tidal forces on the masses (the mirror and the beam-splitter in the case of an 
interferometer, see figure 1). 

A detailed analysis of the gauge of the local observer is given in ref. [23], 
sect. 13.6. Here wc recall only the more important features of this gauge: 

the time coordinate xa is the proper time of the observer O; 

spatial axes arc centred in O; 

in the special case of zero acceleration and zero rotation the spatial coor- 
dinates the proper distances along the axes and the frame of the local 
observer reduces to a local Lorentz frame: in this case the line element reads 

ds^ = -{dx"f + S.jdx'dx^ + 0{\x^'^)dx°'dx'^; (38) 

the effect of SGWs on test masses is described by the equation for geodesic 
deviation in this gauge 
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Photon 




Figure 1: photons can be launched from the beam-sphtter to be bounced back 
by the mirror 



= -m,kox\ (39) 

where Rq^.q are the components of the linearized Riemann tensor [24j . 

Recently, Capozzicllo and Corda [13] have shown that the response of an 
interferometer to SGWs is invariant in three different gauges well known in 
literature and for all frequencies of SGWs, while in previous literature this gauge 
invariance was only proved in the long wavelength approximation [16| . In |13| , 
the angular pattern of interferometers was computed in the TT gauge, while in 
this paper it is computed in the full angular and frequency dependences directly 
in the gauge of the local observer. An analysis similar to the one in refs. [131 [25] 
will be used. We emphasized that such an analyses has been performed for the 
first time in [5S], but it has been strongly generalized to angular dependences 
and massive waves in [13) . 

In the gauge of the local observer, two different effects have to be considered 
in the calculation of the variation of the round-trip time for photons. This works 
in analogy with the cases in refs. [T31[2S] where the effects considered were three, 
but the third effect vanishes putting the origin of the coordinate system in the 
beam splitter of the interferometer. Then equations (74) and (75) in [T3], that 
represent the variations of the coordinates of the mirror of the interferometer in 
presence of a SGW, can be rewritten in the frame of the local observer as 

5x{t) = ixo$(0 (40) 

and 

5y{t) = \x^m- (41) 

Eqs. ([40]) and (|4T]l have been obtained using the perturbation method [24] 
from the equation of motion which are 
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x = ^<i>x (42) 

and 

y = (43) 

For the gauge invariance of the hnearized Riemann tensor [24], eqs. ([42]) and 
can be obtained directly from the TT gauge of eq. ([57)) . 

To compute the response function for an arbitrary propagating direction 
of the SGW one recalls that the arms of the interferometer are in the it and 
if directions, while the x, y, z frame is adapted to the propagating SGW. For 
a better definition of the x^y^z frame we can write, to first order in $, the 
coordinate transformation x" = x°'(x^f) from the TT coordinates (here labelled 
ttt, Xtt,ytt, ztt) to the frame of the local observer as |32j 

t = tu 

X ^ Xtt + ^Xtt'^ 

(44) 

y = ytt + |ytt$ 

Z = Ztt- 

Then, a spatial rotation of the coordinate system has to be performed: 
u = —X cos9 cos(f> + y siiKj) + z siii9 cos(f> 

V = —XCOS0 smcf) ~ y cos(f> + zsiiiO siiKp (45) 

w = xsinfl + ZCOS0, 

or, in terms of the x, y, z frame: 

X ~ —ucosO cos(f> — V COS0 smcf) + w sinO 

y = u sin (j) — v cos 4> (46) 

z = w sin0 cos + w sin0 sin (/) + It; cos 0. 

In this way, the SGW is propagating from an arbitrary direction to the 
interferometer (see figure 2). Assuming that the mirror of eqs. PO)) and (|4ip 
is situated in the u direction, if one uses eqs. (TO)) . (gOl) and ([iT|) . the u 

coordinate of the mirror is 

u = L + ]^L A^{t-us\YL 9 cos (j)). (47) 

where 

A = cos^ cos^ + sin ^ 0. (48) 
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Figure 2: a SGW propagating from an arbitrary direction 



A good way to analyse variations in the proper distance (time) is by means 
of "bouncing photons" (see [13l[25] and figure 1). 

In this Section, a photon which propagates in the u axis is considered, but, in 
next Section, the analysis will be almost the same for a photon which propagates 
in the v axis. 

Putting the origin of the coordinate system in the beam splitter of the inter- 
ferometer and using eq. (j47|) , the unperturbed coordinates for the beam-splitter 
and the mirror are ujy = Q and Um = L. Thus, the unperturbed propagation 
time between the two masses is 

T = L. (49) 

From eq. (|47p the displacements of the two masses under the influence of 
the SGW are 

5ub{t) = (50) 

and 

5u,n{t) ^ ]^LA'^{t - L sine cos (j)). (51) 
In this way, the relative displacement, which is defined by 

SL{t)^Su^{t)-5ubit), (52) 
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gives 



= ^ h A^t-L sine cos (53) 

But, for a large separation between the test masses (in the case of Virgo the 
distance between the beam-splitter and the mirror is three kilometers, four in the 
case of LIGO), the definition ([5^ for relative displacements becomes unphysical 
because the two test masses are taken at the same time and therefore cannot be 
in a casual connection [13l [25] . The correct definitions for the bouncing photon 
are 

5Li{t)=5urrXt)-5ub{t-Ti) (54) 

and 

5L2{t) = 5um{t - T2) - 6ubit), (55) 

where Ti and T2 are the photon propagation times for the forward and return 
trip correspondingly. According to the new definitions, the displacement of one 
test mass is compared with the displacement of the other at a later time to allow 
a finite delay for the light propagation [13l |25] . The propagation times Ti and 
T2 in eqs. ([54]) and (|55|) can be replaced with the nominal value T because the 
test mass displacements are already first order in <I> [13]. Thus, the total change 
in the distance between the beam splitter and the mirror, in one round-trip of 
the photon, is 



6Lr.t.{t) = SLi{t -T)+ 5L2{t) = 25u„,{t - T) - 5ub{t) - Sut{t - 2T), (56) 
and in terms of the amplitude of the SGW: 

SLr.tXt) ^ LA^{t ~ Lsine coscf) - L). (57) 

The change in distance (|57|) leads to changes in the round-trip time for photons 
propagating between the beam-splitter and the mirror: 

= A$(t-Lsin6lcos(/)-L). (58) 



4 Effect of curved spacetime 

In the last calculation (variations in the photon round-trip time which come from 
the motion of the test masses inducted by the SGW), we implicitly assumed that 
the propagation of the photon between the beam-splitter and the mirror of the 
interferometer is uniform as if it were moving in a flat space-time. But, the 
presence of the tidal forces indicates that the space-time is curved (TSj [25] . As 
a result one more effect after the first discussed has to be considered, which 
requires spacial separation [131 125] . 
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From equations ([1^ . (|43p . and the tidal acceleration of a test mass 
caused by the SGW in the u direction is 

u{t — u sin 6* cos (/)) = ^LA^{t — w sin cos (/)) . (59) 



Equivalently one can say that there is a gravitational potential [131 Ell ES] : 

1 r ■■ 

V{u,t)^ — LA I <S>{t- I sine cos 4>)dl, (60) 

2 Jo 

which generates the tidal forces, and that the motion of the test mass is 
governed by the Newtonian equation 

9 = -^V. (61) 

For the second effect, the interval for photons propagating along the u - axis 
can be written like 

ds^ = goodt^ ~ du^. (62) 

The condition for a null trajectory, i.e.ds = 0, gives the coordinate velocity 
of the photons 

v'^{^r = l + 2V{t,u), (63) 
which, to first order in $, is approximated by 

v^±[l + V{t,u)], (64) 

with + and — for the forward and return trip respectively. Knowing the 
coordinate velocity of the photon, the propagation time for its travelling between 
the beam-splitter and the mirror can be defined: 

Ti{t) = / — (65) 

J«,(t-Ti) V 

and 

T.(„=r"' tM. (66) 

The calculations of these integrals would be complicated because the Um 
boundaries of them are changing with time: 

ut,{t) = (67) 

and 



,{t) = L + 6u,n{t). (68) 
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But, to first order in $, tliese contributions can be approximated by 6Li{t) 
and 5^2 (i) (see eqs. ([54)) and ((55|) V Thus, the combined effect of the varying 
boundaries is given by 6iT{t) in cq. ([55]) . Then, only the times for photon 
propagation between the fixed boundaries and L have to be computed. Such 
propagation times will be indicated with ATi 2 to distinguish from Ti_2- In the 
forward trip, the propagation time between the fixed limits is 

r^jH^^ - ['^V{t\u)du, (69) 
Jo v{t',u) Jo 

where t' is the delay time (i.e. t is the time at which the photon arrives in 
the position L, so L — u = t — t') which corresponds to the impcrturbcd photon 
trajectory: 

t' = t- (L-u). 
Similarly, the propagation time in the return trip is 

AT2{t)=L-[ V{t',u)du, (70) 



where now the delay time is given by 

t' ^t-u. 

The sum of ATi{t — T) and A.T2{t) gives the round-trip time for photons 
travelling between the fixed boundaries. Then, the deviation of this round-trip 
time (distance) from its unperturbed value 2T is 



52T{t) = - [Vit -2L + u, u)du+ 
~ llvii-'>J'^iJ')]du, 

and, using eq. ([60)) . 

62T{t) = A J^[J^' 'Pit -2T + l{l-sm0 cos (j>))dl+ 

" lo ^ + sin 6 cos (j))dl]du. 
Thus, the total round-trip proper time in presence of the SGW is: 



(71) 



(72) 



Tt = 2T + 5iT + S2T, (73) 

and 

(5T„ = Tt-2T = SiT + S2T (74) 

is the total variation of the proper time for the round-trip of the photon in 
presence of the SGW in the u direction. 

Using eqs. ([55]). ([7^ and the Fourier transform of <I> defined by 
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/oo 
dt<i>{t)cxp{iujt), (75) 

the quantity (j74p can be computed in the frequency domain as 

6Tu{lo) = SiT{uj) + 62T{lo) (76) 

where 

SiT{u;) = exp[iLuL{l + sine cos (t))]LA^{ui) (77) 



X 'T'/ \ LA[ — l+cxp[iu;L(l+sin 6 cos 0)1 — iLajfl+sin 6 cos 0) . 
d2l (UJ) (l+si„ecos0)2 + 

I cxp(2ia;L) (1— cxp['icjL( — 1+sin ^ cos (p)]-\'iLLU~{l-\-s'm 9 cos 0) i ^ / \ 
+ (_n.sinecos0)2 Jy'l'^j- 



(78) 



In the above computation, the derivative and translation theorems on the 
Fourier transform have been used. In this way, using eq. (|48p . the response 
function of the u arm of the interferometer to the SGW is obtained 

^"(^) = §g - 2i^[-l + exp(2zc.L)+ 

(79) 

+ sin cos (/)((! + exp{2iujL) — 2expiujL{l + sin 9 cos (f>))], 

which is the same result of equation (148) in |13| . where the computation 
has been performed in TT gauge. 

5 Computation for the v arm 

The computation for the v arm is similar to the one above. The coordinate of 
the mirror in the v arm is: 

v^L+^LB<^{t-v sine sm(j>), (80) 

where 

B EE cos^ sin^ (j) + cos^ 0. (81) 

Thus, with an analysis similar to the one of previous Sections, the variation 
in the photon round-trip time which comes from the motion of the beam-splitter 
and the mirror in the v direction is: 

^i^^ = LB$(t-L sin 61 sin 0-L), (82) 
while the second contribution (propagation in a curved spacetimc) is 

S2T{t) = ^LB [J^<i>{t-2T + l{l- sine sin (l)))dl+ 

(83) 

- Jp" $(t - 1{1 + sine sin (f>)dl]du. 



14 



Then, the total response function in the v direction for the SGWs is 



(84) 

+ sin 6' sin (/)((! + exp{2iLLjL) — 2expiujL{l + sin 6* sin 0))], 

which is exactly the result (149) in p3] where the computation has been 
made in the TT gauge. 



6 The total response function 

The total response function is given by the difference of the two response func- 
tions of the two arms: 



Htotioj) = Hu(uj) - H^iuj), (85) 
and, using eqs. ([75]) and ((M|) . one gets 

Htotito) = = i^{cos (/)[l+exp(2iwL)- 2 expitjL(l + sin 61 cos (/))] + 



— sin + exp(2iwL) — 2 exp iujL{l + sin 6 sin (/))]}■ 

(86) 

This equation gives exactly the total response function (150) in [13] where 
the computation has been made in the TT gauge. Thus, the gauge invariance 
between the gauge of the local observer and the TT gauge for the total an- 
gular and frequency dependences of the response function of an interferometer 
has been proved, improving the result in |13| . where the gauge invariance was 
shown only in its simplest geometric configuration (propagation of the SGW 
perpendicular to the two arms). 

Eq. ([86]) is also in perfect agreement with the detector pattern in [l6j and 
|31j in the low frequencies limit: 

Htot (w 0) = - sin^ 0cos2(f>. (87) 



7 The signal to noise ratio in the two antennas 
correlation for the detection of a stochastic 
background of scalar waves 

Let us consider a stochastic background of SGWs (see [lH [22l [30l [33] and [26] - 
[23]). It can be descripted in terms of the scalar field <I> and characterized by a 
dimensionless spectrum (an analogous definition exists for tensorial waves, see 
[26] for example) 
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where 

is the (actual) critical density energy of the Universe [TSl [221 ISHl [S3] and 

I2S]-[2S1. 

Interacting with a stochastic background, the complex Fourier amplitude $ 
can be treated as a random variable with zero mean value in a way similar to 
in the Fourier domain [26l [33] . By assuming that the stochastic background 
is isotropic and stationary, the ensemble average of the product of two Fourier 
amplitudes can be written as 

< $*(/, nw, n') s{f - f')d^Ci, n')s^{f), (90) 

where is a unit vector specifying the propagation direction, and, using the 
explicit definition of the spectrum ([SS]) [26] [33] , 

^'*(/) - • (91) 

The optimal strategy for a potential detection of a stochastic background 
requires the cross correlation of at last two detectors with uncorrelated noises 
ni{t), i ~ 1,2 [26] I33j. Given the two outputs over a total observation time T, 

s,{t) = SUt)+n,it), (92) 
a "signal" S can be constructed: 

S= f ^ si{t)s2{t')Q{t-t'), (93) 

where Q is a suitable filter function, usually chosen to optimize the signal 
to noise ratio (SNR) [M] [33] 

SNR 5 > /A5. (94) 

In the above equation AS* is the variance of S. By assuming that the obser- 
vation time is much larger than the temporal distance between the two detectors 
for which Q is not zero, in the frequency domain it is 

,S>=I^T.Re\[il'S!MB^pmil\ (95, 

where 0{f) is the Heaviside step function and 7(/) the so - called overlap 
reduction function defined in [M] and adapted to scalar waves following [211 [33] ■ 

For the computation of the variance, one assumes that, in each detector, the 
noise is much greater than the strain due to SGWs, obtaining 

A52^-/ d/Pl(|/|)P2(|/|)|Q(/)|^ (96) 



16 



where Pi(|/|) is the one-sided power spectral density of the i detector [251155] . 
Introducing the inner product 

(a, 5) = |^°° d/a(/)fe(/)A(/)P2(/)| , (97) 
the squared SNR can be rewritten as 

^2 (g 



{SNRf = 2T{^Y , A A ■ (98) 



^2 _ c,rr,tlQ ^ \^^ PPl{\f\)P 2{\f\) 

The above ratio is maximal for 



and with this optimal choice the signal to noise ratio becomes 



8 The generalized overlap reduction function for 
the advanced LIGO-LIGO correlation 



The overlap reduction function for SGWs has been used in [TBI ESI [33] starting 
by its definition for tensorial waves in [34| . The expression 

= dClcxp{2mfn ■ r^2)Hi{.f)H-2{.f) (101) 

can be used, where T*i2 is the distance between the two detectors and Hi{f) 
the angular pattern of the i detector (i = 1, 2) [T21[H1[33| In literature, the low 
frequencies approximation angular pattern (|87p has been used in the compu- 
tation of the overlap reduction functions for stochastic backgrounds of SGWs 
[TBI [55] . Now, the analysis will be improved with the aid of the frequency de- 
pendent angular pattern (|86[) . Putting the origin of the coordinate system in 
the LIGO site in Hanford it is 

7(/) = £ / smeded<j)cxp{27rifX)^{cos(j)[l + cxp(2iLjL) - 2 cxp iuL {1 + sin 6 cos cj})]- 

— sm(j)[l + cxp{2iujL) — 2 exp iujL{l + sin 6 sin (/))]}• 

""^^fj~/^^ {cos((/) - + exp{2iujL) - 2cxpiwL(l + sin(6l - O^) cos{(j) - 0i))] + 

— sin(0 — (/)i)[l + cxp(2iLi;L) — 2expiwL(l + sin(6' — 6i) sin(0 — 0i))]} 

(102) 

where 0i = -28M^,(f>i = 71.2", L = AKm and 
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Figure 3: The absolute value of the overlap reduction function (|102p in the 
frequency-range IQHz < / < IKHz. 



X = (i{cos(0i — sin6'i + cos 6*1 (cos + sin(/)i))- 
(— sin 6* + cos^(cos + sin (/))) + (cos^i + sin (cos0i + sin 0i)) • (103) 
(cos0 + sn\0{cos(f) + sin(/))) — sin(0i + (f) 

with d = 2997. 9A'm. For the position and orientation of the two LIGO sites 
ref. [36] has been used. 

In figure 3 the absolute value of the overlap reduction function (jl02p is drawn 
in the frequency-range of Earth based interferometers that is the interval 

IQHz < f < IKHz. 

In the figure the value of the overlap reduction function of the two LIGO 
interferometers for SGWs appears low, thus, in principle, we need a long inte- 
gration time to improve the SNR. A lower bound for the integration time of 
a potential detection will be computed in next Section. This lower bound for 
the integration time corresponds to an upper bound for the specific SNR of the 
cross-correlation between the two advanced LIGO. 
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9 The specific signal to noise ratio for scalar 
waves in the two advanced LI GO interferom- 
eters 



For the computation of the SNR (|100p the analytical fit of [37] for the noise 
spectral density of advanced LIGO is used: 

pif) = mi)-'-'' ~ Hir' + '''^\fSl;^^^^^ ] f>fs 

P{f) - oo / < /, 

where Pq = 6 • 10-^^Hz~^, fo = 215Hz and /, = 20Hz. 

The spectrum of stochastic backgrounds is flat in the frequency-range of 
Earth based interferometers while it is well known that WMAP observations put 
strongly severe restrictions on the same spectrum HU [3D1 1311 |3H] and 
|29| . In fig. 4 we map the spectrum Qsgw choosing the amplitude (determined by 
the ratio j^^^ [20l[2Tl[30l[35l [38]) to be as large as possible, consistent with the 
WMAP constraints [20l [211 [30l [35l [38] . The inflationary spectrum rises quickly at 
low frequencies (wave which re-entered in the Hubble sphere after the Universe 
became matter dominated) and falls off above the (appropriately redshifted) 
frequency scale fmax associated with the fastest characteristic time of the phase 
transition at the end of inflation. The amplitude of the flat region depends on 
the energy during the inflationary stage [|ol [HI [SO] [35l [38] . Because WMAP 
data are consistent with a maximum inflationary scale Mini = 2 • lO^^GeV 
[101 [m [301 [33 [Ml, this means that today, at LIGO and LISA frequencies, 
indicate by the arrows in fig. 4, 



^^s3«,(/)/^?oo <9*10-". (105) 

Using equations ()100|) . (|102p . (|104|) and ()105p the cross-correlation needs a 
time of ~ 8 * lO^years in order to obtain (SNR) = 1. 

We emphasize that the assumption that all the scalar perturbation in the 
Universe are due to a stochastic background of SGWs is quit strong, but the 
result can be considered like a lower bound for the integration time. 

To better understand the difference between the response function of this 
paper and previous low-frequency approximated one a brief comparison is now 
performed. 

In the computation of the time that the cross-correlation needs to obtain 
{SNR) ~ 1 one can replace the response function ([86]) with the approximate 
one l[87]) in cq. (jlOOp . Integrating in all the frequency range of Earth based 
interferometers, the underestimated time ~ 3 * 10^ years is obtained (see [20]). 
This is because the correct response function ([55]) drives the overlap reduction 
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function to a stronger decrease with respect to the constant value of the ap- 
proximate response function ([87|) . The correct response function ([86)) in the 
Hanford - Livingston direction is shown in figure 5. One could also replace the 
response function ([5S)) with the approximate one ([57)) in eq. (jlOip . integrating 
only in the low frequencies range 10 Hz < / < lOOHz, where the approximation 
is correct. In that case the result is overpriced: ~ 7 * lO^years. 

However the three results (the correct one and both the underestimated and 
the overpriced ones) are too much long. This means that for a detection we 
hope in the LISA interferometer and in a further growth in the sensitivity of 
advanced projects. In this second case one notes from eq. (jlOOp that only a 
gain of 5 orders of magnitude in the noise spectral density of advanced LIGO 
could in principle provide (SNR) ~ 1 for integration times which are not too 
long (order of some years). 



20 



Ene 



rgy 



LISA 



VIRGO 



6 



Hz { - Graphics - ) 



15 



10 



5 



5 



10 




10 



12 




14 



-15 - 



18 



Figure 4: The spectrum of relic SGWs in inflationary models is flat over a wide 
range of frequencies. The horizontal axis is logj^Q of frequency, in Hz. The verti- 
cal axis is logj^o ^gsw- The inflationary spectrum rises quickly at low frequencies 
(wave which rentered in the Hubble sphere after the Universe became matter 
dominated) and falls off above the (appropriately redshifted) frequency scale 
fmax associated with the fastest characteristic time of the phase transition at 
the end of inflation. The amplitude of the flat region depends on the energy 
during the inflationary stage; we have chosen the largest amplitude consistent 
with the WMAP constrains on scalar perturbations. This means that at LIGO 
and LISA frequencies, ilsgw{f)hiQQ < 9 * 10"^'^. Adapted from C. Corda - Pri- 
mordial production of massive relic gravitational waves from a weak modification 
of general relativity - Astropart. Phys. 30, 4, 209-215 (2008) 



10 Conclusions 

With the "bouncing photon" analysis, in [13] the gauge invariancc of the re- 
sponse of an interferometer to SGWs has been recently demonstrated in three 
different gauges well known in literature and for all the frequencies of SGWs, 
while in previous literature this gauge invariance was only shown in the long 
wavelength approximation. In this paper the analysis of the response function 
for SGWs has been generalized in the full angular and frequency dependences 
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Figure 5: the absolute value of the total response function of the LIGO inter- 
ferometer sto SGWs for the Hanford - Livingston direction. 
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and directly in the gauge of the local observer, which is the gauge of a laboratory 
environment on Earth. The result has been used to analyse the cross-correlation 
between the two LIGO interferometers in their advanced configuration for a po- 
tential detection of a stochastic background of SGWs and to release a lower 
bound for the integration time of a potential detection. A comparison between 
the new response function and previous low-frequency approximated ones has 
also been performed. Because the lower bounds result very long, for a detection 
we hope in the LISA interferometer and in a further growth in the sensitivity 
of advanced projects. 
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